Astronomy & Astrophysics manuscript no. 13885TEX 


© ESO 2010 


March 30, 2010 





Properties and nature of Be stars * 

27. Orbital and recent long-term variations of the Pleiades Be star 

Pleione = BU Tauri 

J. Nemravova 1 , P. Harmanec 1 , J. Kubat 2 , P. Koubsky 2 , L. Iliev 3 , S. Yang 4 , J. Ribeiro 5 , M. Slechta 2 , L. Kotkova 2 , 

M. Wolf 1 , and P. Skoda 2 



o 

(N 



(N 

Oh 
6 



> 

in 

in 

co 
o 
o 



X 



Astronomical Institute of the Charles University, Faculty of Mathematics and Physics, 
V Holesovickach 2, CZ- 1 80 00 Praha 8, Czech Republic 

Astronomical Institute of the Academy of Sciences, CZ-251 65 Ondfejov, Czech Republic 

Institute of Astronomy, Bulgarian Academy of Sciences, BG-1784, 72 Tsarigradsko Chaussee Blvd., Sofia, Bulgaria 
Physics & Astronomy Department, University of Victoria, PO Box 3055 STN CSC, Victoria, BC, V8W 3P6, Canada 
Observatorio do Instituto Geografico do Exercito, R. Venezuela 29, 3 Esq. 1500-618, Lisboa, Portugal 



Release March 30, 2010 



ABSTRACT 



Radial-velocity variations of the Ho- emission measured on the steep wings of the Ho- line, prewhitened for the long-time changes, 
vary periodically with a period of 218 c !025±O!022, confirming the suspected binary nature of the bright Be star BU Tau, a member 
of the Pleiades cluster. The orbit seems to have a high eccentricity over 0.7, but we also briefly discuss the possibility that the true 
orbit is circular and that the eccentricity is spurious owing to the phase-dependent effects of the circumstellar matter. The projected 
angular separation of the spectroscopic orbit is large enough to allow the detection of the binary with large optical interferometers, 
provided the magnitude difference primary - secondary is not too large. Since our data cover the onset of a new shell phase up to 
development of a metallic shell spectrum, we also briefly discuss the recent long-term changes. We confirm the formation of a new 
envelope, coexisting with the previous one, at the onset of the new shell phase. We find that the full width at half maximum of the 
Ha profile has been decreasing with time for both envelopes. In this connection, we briefly discuss Hirata's hypothesis of precessing 
gaseous disk and possible alternative scenarios of the observed long-term changes. 

Key words, stars: early-type - stars: binaries - stars: Be - stars: individual: BU Tau 



1. Introduction 

Pleione (BU Tau, 28 Tau, HD 23862) is a well-known Be star 
and a member of the Pleiades cluster. It underwent several 
phase transitions between B, Be, and Be shell phases, accom - 
panied by prono u nced li g ht variation s; see, e .g. Gulliverl(fT9T7l). 
ISharov & Lvutvl (Il976l). IllievetalJ dl988h. Isiiarov & Lvutvil 
d 1992b iHirata & Kogurel (1 1976b . IHirata & Ko surd (119771) . 
Hiratal d 1995b. iDoazan et all d 1988b . IllievetalJ (120071) . and 



Tanaka et al 



(120071) . 
rather 



There is a rather complicated history of attempts to 
study the radial-velocity (RV hereafter) variations of this star. 
IStruve & Swingsl dl943i) measured RVs on the photographic 
spectra taken in the years 1938-1943 and tentatively concluded 
that the RV of BU Tau va ries with a pos sible period of 142 days 
or- less likely - 106 davs. lMerrillldl952b studied RVs from 1941 
to 1951 and found clear long-term vari ations wi t h som e overlap- 
ping changes on a shorter time scale. iGulliverl (1 1977b analyzed 
a large collection of digitized photographic spectra from 1938- 
1954 and from 1969-1975 and concluded that there are no sig- 
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nificant RV changes. iBallereau et al.l (fl988) carried out an anal- 
ysis of a homogeneous series of Haute Provence high-dispersion 
photographic spectra from 1978-1987 and once more concluded 
that the shell RVs var y with periods of 136.0 and 106.7 days. 
Kata hira et al .1 ( 1996a b) analyzed shell RVs from the two con- 
secutive shell phases separated some 34 years, using published 
as well as new RVs and concluded that BU Tau is a spectro- 
scopic binary with an orbital period of 218^0, semi-amplitude of 
5.9 km s" 1 , and a la rge orbital eccentricity of 0.60. However, 
iRivinius et al.l (12006b - analyzing a series of electronic spec- 
tra - were unable to confirm the 218-d pe riod and conc luded 
that BU Tau is not a spectroscopic binary. iHiratal ^2007) ana- 
lyzed a long series of polarimetric observations and presented 
a model of a slowly precessing disk to explain the long-term B 
- Be - Be shell phase transition. He argued that the disk pre- 
cession is caus ed by the attrac tive force of the secondary in the 
218-d binary. Harmanec ( 1982) compiled the majority of at that 
time available RVs of BU Tau and averaged them over about 
100 days. This resulted in a smooth RV curve with a period 
of about 1300 days (35.6 year s), in phase with the recorded 
shell episodes. lHarmane"cl d!982h speculated that BU Tau could 
be a long-periodic binary with shell phases occurring always at 
the same orbital phases. A more distant companion with an an- 
gular distance of 0722 was indeed discovered from speck le in- 
terferometry bv lMcAlister et all (Il989b . iGies etalJ (1 1990b stud- 
ied a sequence of low-dispersion Ha spectra of BU Tau taken 
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Fig. 1. Comparison of Ha profiles from different stages of the 
long-term changes. 



with a sampling rate of 7 ms during a lunar occultation on 
1987 March 6. They detected an asymmetry of the envelope in 
agreement with the observed long-term VI R changes. They spec- 
ulated that the speckle-interferometric component could have 
an eccentric orbit and that the rec urrent shell phases could 
be ca used by its periastron passages. iLuthardt & Menchenkoval 
( 1994) compiled RVs from the years 1938-1990 and confirmed 
a period of 12450-12860 days. They advocated an eccentric or- 
bit and mass transfer resulting in a release of a new shell during 
periastron passages, but the gaps in their RV curve do not allow 
one to conclude that the orbit has a high eccentricity. Finally, 
usin g the technique of ad aptive optics photometry and astrom- 
etrv. iRoberts et al.l d2007l) report discovery of a new companion 
to BU Tau at a separation of 4'.'66 with a spectral type M5. They 
also confirm a companion at 0724 and discuss other suggested 
companions. 



Table 1. Journal of new spectroscopic observations for BU Tau. 



Station 


Time interval 


No. of. 


Wavelength 


Source 


(HJD-2400000) 


obs. 


region (A) 


Ondfejov 


49581 -54872 


101 


6200 - 6800 


DAO 


49786 - 54912 


26 


6150 - 6700 


OHP 


51569 - 52664 


21 


6200 - 6700 


Rozhen 


52710 - 54108 


23 


6520 - 6610 


Lisboa 


54874 - 54881 


4 


6520 - 6600 




Fig. 2. A complete series of our Ha profiles in a gray repre- 
sentation (only a few saturated or underexposed spectra were 
omitted). Abscissa shows the wavelength scale in A, while the 
time on ordinate is shown in JD-2400000. Each horizontal strip 
represents an average of spectra secured within 200 days, and 
dark horizontal belts correspond to time intervals from which no 
spectra are available. At the bottom, there is a scale showing the 
correspondence between the flux level in the units of continuum 
and the gray scale. 



1 1 1 1 


v f HJ] 






. f H. 


D=245471 8.51 38 




' x -^~ H 


-24541 1 5 l.lSfi 












= 2453335.5298 




~ — • — — '~k7 


= 24b2877.5722 



6674 6676 6678 6680 668? 6684 6686 6688 6690 
wavelength(A) 



We succeeded in collecting a rich series of electronic spectra 
at several observatories, covering many cycles of the suspected 
218-d period. The main goal of this study is, therefore, to re- 
solve the issu e of wh ether BU Tau is a spectroscopic binary. 
Kata hira et al.l (1 1996 a b) based their orbit on the RV measure- 
ments of shell lines that may be affected by possible asymme- 
tries in the circumstellar matter. Moreover, their RV curve has 
a rather small amplitude and is based on a collection of hetero- 
geneous data. It naturally shows a rather large scatter around 
the mean curve. The spectra at our disposal all cover the red 
spectral region near Ha. They were taken over the time inter- 
val when the star had fairly strong Ha emission. Therefore, 
our study is based on the RV measurements of the steep wings 
of the emission, which is a procedure that turned out to be 



Fig. 3. Selected He I 6678 A line profiles, ordered in time, with 
corresponding HJDs. 



successful for detecting the duplicity of several other Be stars 
dBozic et al.lll995b iKoubsky et al Il2000k lHarmanec et al J 120001: 
iMiroshnichenko et al.ll2001ll2002l) . 

Since very pronounced long-term spectral variations oc- 
curred over the time interval covered by our spectra, we also 
briefly describe these changes and discuss them, especially in 
relation to the model put forward by Hirata (20071). 
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HJD - 2400000 

Fig. 4. Measured RVs of the Ha emission wings plotted vs. time. 
Prewhitening for the long-term changes, carried out with the 
help of the program HEC13, is shown by a line. Empty squares 
show the alternate way to remove long-term RV changes via in- 
dividual y velocities for subsets spanning no more than a year. 
See the text for details. 
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Fig. 6. A comparison of two pairs of the Ha line profiles from 
the locally recorded velocity extrema (HJDs of the profiles are 
indicated). 
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Fig. 5. Measured Ha absorption core plotted vs. time. We also 
include d the RV measurements o f shell lines by Katahir a et al.l 
(Il996bl) and lRivinius et al.l (12006) to this plot. Prewhitening for 
the long-term changes, carried out with the help of the program 
HEC13, is shown by a line. See the text for details. 

2. Spectroscopic observations and their 
reductions 

The red spectra at our disposal were obtained at five observato- 
ries and their overview is in Table Q] Details about the instru- 
ments and data reduction can be found in Appendix [A] where 
also Table I A. ll with our RV measurements of the steep wings of 
the Ha emission and of the Ha absorption core is provided. The 
latter was m easured for compariso n with the RVs collected and 
analyzed by Katah ira et al.l d!996bl) . but only for those spectra 
where the absorption was clearly visible. 

Over the interval of the more than 5000 days covered by our 
observations, the strength of the Ha emission gradually declined 
and the shape of the Ha profile underwent notable changes. 
Typical examples for several distinct stages are shown in Fig.Q] 
and the whole development of a new shell and metallic-shell 
phase is shown as a gray-scale representation of all usable Ha 
profiles in Fig. [2] 



The fading of the Ha emission was accompanied by a light 
decrease in the J, H, K, and L IR photomet ric bands that started 
around JD 2451500 dTaranova et al J 12008). This clearly corre- 
sponds to the gradual development of the hydrogen she l l spec - 
trum according to our spectra - see also iTanaka et al.l d2007l) . 
Emission has been slowly fading from JD 2453000 until now, 
when its peak intensity represents only about 30% of the in- 
tensity seen in our earliest spectra. During the transition from a 
single-peaked to double-peaked emission, there is some time in- 
terval when the Ha profile has a characteristic wine-bottle shape. 
The occasional presence of additional absorption components 
has be en already noted by llliev et ail (120071) or ITanaka et al.l 
(200?]) and is typical of all recorded shell phases of BU Tau. 
Besides the occasional presence of one or more additional ab- 
sorptions, extended red emission wings are seen on some Ha 
profiles. This makes the emission wings asymmetric and hard 
to measure for RV. We also note that all double-peaked profiles 
recorded prior to about JD 24540000 always have a red peak 
stronger than the violet one. Figure [2] shows that the width of the 
Ha emission has remained more or less constant over the whole 
time interval covered by our observations. The same figure also 
shows that the metallic shell phase appeared rather abruptly. 

Figure [3] shows the gradual development of the He I 6678 A 
line profile. It illustrates well how shallow the line is at the be- 
ginning of a new shell phase. A very interesting finding is that, 
even for the B8 star, a presumably photospheric He I line can 
develop a shell component. The profile clearly gets stronger and 
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Fig. 7. Stellin gwerl dl978l) PDM 6 statistics for all emission- 
wing _RVs__(jto2) a nd she ll absorption-core RVs includ- 
ing iKatahira et al.l (fl996b). The dominant frequency of 
0.004587 c d -1 corresponds to the 218-d period. 



narrower as the hydrogen shell line gets deeper. The additional 
absorption at the blue wing of the line seen on more recent spec- 
tra is the Fe II 6677.305 A shell line. 



3. Radial-velocity changes 

Figures |4] and [5] are the time plots of the measured RVs vs. 
time for the Ha emission wings and the absorption core. In 
the later, we a l so inclu ded all shell RVs u s ed an d published by 
IKatahira etail d!996bl) and iRivinius et all d2006l) . One can see 
systematic RV changes on at least two distinct time scales: a 
smooth change on a longer time scale and overlapping more 
rapid changes, especially the occasional steep decreases in RV. 

Considering the uncertainties in accurate RV measurements 
combined with the fact that the full amplitude of the changes 
is low, it was deemed useful to convince readers that the RV 
changes are not only a result of changing asymmetry of the pro- 
files, but they also represent a real shift of the whole line. To this 
end, we compare in Fig. |6]two pairs of the Ha line profiles ob- 
tained near the local RV extrema. The upper pair comes from the 
beginning of a new shell phase and the bottom one from a more 
recent time when a weaker emission and deeper shell cores are 
present in the profiles (note a large difference in the flux scale of 
the two plots). The RV shift of the whole emission and absorp- 
tion core is seen beyond any doubt. We, therefore, conclude that 
our RV measurements reflect real RV variations of BU Tau. 



In accordance with Katahira et al.l (Il996bl) . we find that the 
evolution of the emission episode is accompanied by long-term 
RV changes that need to be removed prior to a search for possi- 
ble periodic RV changes. To also make this step as objective as 
possible, we used two different procedures. 

One is that we smoothed the long-term changes using the 
program HEC13, w ritten by PH and based on a smoothing tech- 
nique developed by Vondrak d 19691 [l977)0] For both emission 
and absorption RVs, optimal smoothings were obtained for the 
smoothing parameter s = 10~ 16 fitted through 200-d normals. 
(Inspecting the time plots of RVs, we identified ~ 200 days as 
a time scale on which more rapid changes were observed, and 
this was the reason for the choice of 200-d normals. We have 
verified, however, that the result of smoothing is not sensitive to 
the particular choice of the averaging interval for the smoothing 
within reasonable limits.) The RV residuals from the smooth- 
ing we re subjected to a period search based on the Stellingwerf 
(1978) PDM technique over a period range from 5000 down 
to 0.05 d. The dominant frequency found in both searches was 
0.004587 c d~' and its integer submultiples. The one-day aliases 
were largely supressed thanks to having data from observatories, 
that have a large difference in their local time, producing much 
shallower minima in the 6 statistics (~ 0.75 - 0.82) and scat- 
tered phase diagrams. To make the diagrams readable, we show 
the corresponding 6 statistics in Fig.|7]for the emission (top) and 
absorption (bottom) RVs only for a limited frequency interval 
down to 0. 1 c d _1 . The result seems to con firm the 218-d period- 
icity discovered by Kat ahira et al.l d!996bl) . 

As another demonstration that the 218-d period is real, 
we show phase plots in Fig. [8] for the original RVs (without 
prewhitening for the long-term changes) for several subsets of 
data covering time intervals no longer than one year. Clearly 
similar RV curves, with sharp minima, rather flat maxima, and 
a mutual phase coherence, are seen in all cases. The first sub- 
set is based solely on the RVs from the Ondfejov spectra se- 
cured with Re t icon d etector, which were already investigated by 
IRivinius et ail d2006l) . 



4. BU Tau as a spectroscopic binary 

Our findings, and especially the fact that the Ha Balmer emis- 
sion line moves in RV as a whole in spite of very large secular 
changes of its strength, indicate that BU Tau is indeed a single- 
line spectroscopic binary that moves in a highly eccentric orbit. 
We therefore used the program SPEL (written by the late Dr. Jin 
Horn and never published) to derive th e orbital elements. For 
comparison with Katahira et al. ( 1996b), we first derived orbital 
elements for Balmer abs orption RVs, using th e data from their 
study, RVs published by Rivinius et al.l d2006l) . and our own Ha 
absorption RVs, prewhitened with HEC13 as shown in Fig. [5] 
The resulting orbital elements are given as solution 1 in Table [2] 
and the corresponding phase plots are shown in Fig. [9] For more 
clarity, we plot ther e the photographic RVs, Heros RVs from 
Rivini us et al ] (120061) . and our Ha absorption RVs in three sep- 
arate panels. Although R ivinius et al. I d2006h write that the sus- 
pected binary nature of BU Tau could not be confirmed on the 
basis of their data, their RVs also nicely follow the 218-d period. 
This constitutes yet another support for the rea lity of this pe- 
riod. Ou r solution 1 agrees well with the result of lKatahira et al.l 
(fT996bh . 



1 The program HEC13 with brief instructions how to use it is avail- 
able to interested users at http://astro.troja.mff.cuni.cz/ftp/hec/HEC13 . 
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Next, we analyzed the emission RVs that we consider as 
most realistically describing the true orbital motion. To see how 
sensitive the result is to the manner of prewhitening the data we 
derived the elements not only for the RVs prewhitened with the 
help of HEC13 (see above) but also from the original data. To 
this end, we divided the data into subsets spanning no more than 
one year and allowed SPEL to derive separate y velocities for 
individual data subsets. The results are summarized in Tabled 
and the corresponding RV curves compared in Fig.fTUl 

Table 2. Several sets of orbital elements: Solution 1... Katahira, 
Rivinius and this paper, prewhitened with HEC13; Solution 2... 
New emission-line RVs prewhitened with HEC13; Solution 3... 
New emission-line RVs with allowance for locally derived y ve- 
locities. 



Solution: 


1 


2 


3 


Orbital 


Old & new 


Ha emis. 


Ha emis. 


element 


Ha abs. 


wings 


wings 


P(d) 


218.023±0.023 


218.099±0.050 


218.053±0.053 


^"periastr. (m) 


40040.4±1.6 


52039.34±0.69 


52039.73±0.73 


^super.c. (d) 


40032.3 


52035.50 


52034.79 


Tmin.RV (A) 


40044.5 


52040.64 


52041.11 


e 


0.596±0.035 


0.774±0.028 


0.745±0.026 


con 


147.7±4.5 


154.2±4.0 


157.3+3.5 


Ki (km s _1 ) 


5.41+0.35 


6.30±0.63 


6.39±0.46 


y (km s -1 ) 


-0.15+0.1 


0.35±0.15 




rms (km s" 1 ) 


3.21 


1.93 


1.58 



The inspection of Fig. [TUl shows that even the Ha emission- 
wing RVs are indeed indicative of an orbit with high eccentric- 
ity but that there is also an alternative possibility that the ob- 
served deep RV minimum could be a consequence of some un- 
specified effect of circumstellar matter, reminiscent of "an in- 
verse rotational or Rossiter effect". In this case, the true orbit 
could essentially be circular. To this end, we derived yet another, 
a circular-orbit solution for the Ha emission RVs prewhitened 
for long-term changes via HEC13, omitting all RVs from the 
phase interval around phase zero with the most negative RVs. 
This resulted in the following elements: P = 218^34 + 0.62, 
r sup er.c. = HJD 2452009.9 + 4.8, K x = 1.72 + 0.21 km s _1 . 

Using the eliptical-orbit elements for the Ha emission RVs 
from Table [2] we estimated the basic properties of the binary 
from the mass function f(m) = 0.00165 M Q for several plausi- 
ble orbital inclinations, assuming a nor mal mass o f the pr imary 
corresponding to its spectral type after Harmane to be 

Mi = 2.9 M . 

The results of Table [3] show that the binary properties, es- 
pecially the low mass ratio, are quite similar to other binaries 
discovered so far with Be primaries. For the estimates, we only 
considered higher orbital inclinations since BU Tau is one of 
the cases of an inverse correlation between the brightness and 
emission-line strength, whic h indicates that w e see the system 
roughly equator-on - cf, e.g., Har manea(ll983l) . 

If we adopt th e dista nce to Pleiades d - 138 pc after 
iGroenewegen et all d2007l) . we estimate that the projected an- 
gular distance of the binary components should be — 070075, 
dropping down to 070018 at periastron. This angular separation 
is certainly within reach of existing large optical interferometers. 
The only problem is the luminosity ratio primary/secondary. If 
the secondary would be a normal late M dwarf corresponding to 



Table 3. Basic physical properties of BU Tau as a single-line 
binary based on elliptical-orbit solution for the Ha emission RVs 
- cf. Table [2] The estimates are derived assuming the primary 
mass of Mi = 2.9 M , A and A per ; denote the semi-major axis 
and the binary separation at periastron, respectively. 



i 


M 1 /M l 


M 2 


A 


^pen. 


C) 




(Mb) 


(Rg) 


(Re) 


90 


0.0876 


0.254 


223.5 


53.0 


70 


0.0936 


0.272 


223.9 


53.1 


50 


0.1164 


0.338 


225.5 


53.4 



its mass, it would be fainter in the visual region by more than 
10 magnitudes and the only chance to search for it would be in 
the far IR region, where, however, the IR excess from the Be en- 
velope can complicate the detection. However - if it were a hot 
subd warf, similar to the one found for another Be binary <p Per 
bv lGies eta il (119981) - it might be observable in the optical re- 
gion since the absolute visual magnitude of BU Tau is fainter for 
some 2 magnitudes than for the tp Per B0.5e primary. Finally, 
a cool Roche-lobe filling secondary seems improbable since it 
would probably produce binary eclipses. 

In any case, attempts to resolve the 218-d binary system with 
some large interferometer are very desirable since a visual or- 
bit would help not only to estimate the true orbital inclination 
but also to clarify whether the orbit has a high eccentricity or is 
nearly circular. 

5. Comments on Hirata's model 

We have postponed a detailed study of the long-term changes 
for a later work (Iliev et al. in prep.), but we wi sh to commen t 
briefly on the hypothesis put forward recently bv lHiratal d2007l) . 
He obtained systematic spectroscopy and polarimetry of BU Tau 
from 1974 to 2003 and finds a change in the polarization angle 
from about 60° to 130° over that time interval. He interprets this 
change as evidence of the precession of the circumstellar disk 
that is responsible for the observed Ha emission. He further ar- 
gues that also the change in the Ha profiles from a weak double 
emission with a strong central absorption core to a strong emis- 
sion with a wine-bottle shape indicates that the disk was first 
seen m ore or less edge-on and later more face-on. iTanaka et al.l 
(2007) studied the spectra of BU Tau from Nov. 2005 until April 
2007, which cover the period of a formation of the new shell 
phase. They argue that a new disk was formed in the equato- 
rial plane of the B star while the old disk was decaying but 
still present. According to their interp retatio n , the o ld disk was 
precessing in space as suggested by iHiratal d2007l) . Our spec- 
tra co ver a much longer time interval, including the one stud- 
ied by ITanaka et al.l d2007l) . and as Fig. [2] shows, the change of 
the Ha profile was smooth. We thus measured the full width at 
half maximum (FWHM) of a representative selection of our Ha 
emission-line profiles and the variation in FWH M with t ime is 
shown in Fig.Qj] It was already demonstrated bv lStruvel d!931l) 
in his first model of Be stars as rapidly rotating objects that there 
is a clear correlation between the width of presumably photo- 
spheric He I lines and the width of the Balmer emission lines, 
which is preserved during the long-term changes. This correla- 
tion has b een confirmed b y a number of later studies - see, e.g., 
Fig. 5 of ISlettebakl d 1979b . One would therefore expect that, if 
the appearance of a new shell phase of BU Tau is primarily a con- 
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sequence of a geometrical effect, namely a gradual precession of 
a flat disk that becomes to be seen equator-on, the FWHM should 
gradually grow as the new shell phase is approaching. In con- 
trast, Fig. [TT] shows that the FWHM of Ha was slowly decreas- 
ing during the last 15 years. Its dramatic increase is related to the 
formation of a new envelope, which our spectra clearly confirm 
- see Fig. [T2j The apparent discontinous increase in the FWHM 
occurs at the moment when the strength of the broader emission 
from the new envelope rises to a half of the peak intensity of 
the original emission. All this indicates that the observed varia- 
tions are primarily due to physical changes in the circumstellar 
matter and cannot be reduced to a simple geometrical cause - a 
precession of the original gaseous disk. There has been a rather 
widespread tendency in recent years to intepret the presence of 
shell absorption lines as evidence of an equator-on view, since 
many investigators are picturing the Be star disk as a flat struc- 
ture lo cated at the st e llar equator with a (rather sm a ll) opening 
angle (IWatersI 1 19861: iBiorkman & Cassinellil 11993b iHanuschikl 
I1995L 119961) . It is true that this model can lead to theoretical 
Balmer profiles similar to the ob served ones, s ee, e.g., the 3D 
radiative line transfer models bv lHummel {1994). One should be 
aware, however, that there is no unique proof of a specific ge- 
ometry on the level of various simplifications of current models. 
For instance, Hoflich (1987,1988) succeeded in modeling sev- 
eral Balmer emission-line profiles of particular Be stars with his 
model consisting of an NLTE atmosphere and a spherical en- 
velope. It is then conceivable that strong shell lines could also 
develop in the spectrum of a Be star seen more or less pole-on 
in situations where a very extended spheroidal envelope forms 
around it. Similarly, it might be worth considering whether the 
asymmetry detected by the gradual change in the polarimetric 
angle is indeed caused by the precession of a flat disk or by 
some other effect, e.g. by a slowly revolving elongated (non- 
axisymmetric) disk. 
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Fig. 8. Orbital RV curves of the Ha emission shown for subsets 
of data spanning less than a year. For all plots, period 218 c ?053 
was used, with phase zero at HJD 2452041.11, which corre- 
sponds to the RV minimum (see Tabled. 
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Fig. 9. Top: The phase plots of all available Balmer absorption 
RVs, prewhitened for the long-term RV variations with HEC13 
(as shown in Fig. [5j>. Elements from solution 1 of Table [2] were 
used, with phase zero at minimum RV. For clarity, we show three 
diffe rent data subsets separ ately: Top panel: Photographic RVs 
from lKatahira et all dl996b|); Central panel: RVs from electronic 
Heros spectra published bv lRivinius et al.l (12006): Bottom panel: 
RVs from electronic spectra used in this paper. 
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Fig. 11. A time development of the FWHM (in A) of the Ha 
emission. The rapid increase is caused by the formation and a 
fast strengthening of another double emission due to a newly 
formed envelope. 
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Fig. 12. A series of the Ho- profiles over the time interval of the 
formation of a new shell. The HJDs-2400000 of individual spec- 
tra are shown and the time runs from the top to the bottom. One 
can see how the new broad emission gradually rises in inten- 
sity and how its blending with the decaying previous double (but 
narrower) emission creates a profile with four emission peaks for 
some time. Then the new emission gets so strong that it merges 
with the original one. 
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Fig. 10. The orbital RV curves of BU Tau based on the Ha emis- 
sion RVs plotted for the solutions 2 and 3 of Table[2] Phase zero 
corresponds to the respective epoch of minimum RV and the O- 
C deviations from the solutions are shown by small circles in 
separate panels. Top two panels: RVs prewhitened via HEC13 
(solution 2); Two bottom panels: original RVs minus locally de- 
rived systemic y RVs (solution 3). See the text for details. 
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Appendix A: Overview of available spectroscopic 
observations 

Here, we provide some details on the spectra used in this study 
and listed in Table Q] and on their reduction: 

1. Ondfejov spectra: All 101 electronic spectrograms were ob- 
tained in the coude focus of the 2.0-m reflector and have a 
linear dispersion of 17.2 A mm -1 and a 2-pixel resolution 
12600 (11-12 km s _1 per pixel). The first 35 spectra were 
taken with a Reticon 1872RF linear detector and cover a 
spectral region from 6300 to 6730 A. Complete reductions of 
these spectrograms were carried out by JN with the program 
SPEFO, written by the late Dr. J. Horn and further developed 
by Dr. R Skoda and more recently by Mr. J. Krpata - see 
iHorn et al .] (Il996h and lSkodal dl996). The remaining spectra 
were secured with an SITe-5 800 x 2000 CCD detector and 
cover a slightly longer wavelength interval 6260-6760 A. 
Their initial reductions (bias subtraction, flatfielding, cre- 
ation of 1-D images, and wavelength calibration) were car- 
ried out by MS in IRAF. 

2. DAO spectra: These spectrograms were obtained in the 
coude focus of the 1.22-m reflector of the Dominion 
Astrophysical Observatory by SY, who also carried out their 
initial reductions (bias subtraction, flatfielding, and creation 
of 1-D images). Their wavelength calibration was carried 
out by JN in SPEFO. The spectra were obtained with the 
32121H spectrograph with the IS32R image slicer. The de- 
tectors were UBC-1 4096x200 CCD for data before May 
2005 and SITe-4 4096x2048 CCD for data after May 2005. 
They cover a wavelength region from 6150 to 6750 A, have 
a linear dispersion of 10 A mm -1 and 2-pixel resolution of 
21700 ( ~ 7 km s _I per pixel). 

3. OHP spectra: The public ELODIE archive of the Haute 
Provence Observatory (Mou ltaka et al.l 120041) contains 30 
spectra listed as BU Tau, but some of them are actually spec- 
tra of 27 Tau. We were able to recover 21 usable spectra. For 
the purpose of this study, we extracted, rectified, and mea- 
sured only the red parts of these spectrograms. 

4. Rozhen spectra: All 23 spectra from Rozhen observatory 
were obtained in the coude spectrograph of the 2-m RCC 
telescope. A CCD camera Photometries AT200 with SITe 
SI003AB 1024x1024 chip was used. The spectrograph was 
used in a configuration providing high-resolution spectra 
suitable for revealing fine details and the structure of the 
spectral lines. A Bausch&Lomb 632/22.3 grating was used 
in its 2nd order, giving a linear dispersion of 4.2 A/mm 
with 2-pixel resolution of 33000 (~ 4.5 km/s per pixel). 
Wavelength coverage is about 100 A around Ha. The ini- 
tial reduction (bias subtraction, flatfielding, creation of 1-D 
images and wavelength calibration) was carried out by LI in 
MIDAS. 

5. Lisboa spectra: These 4 CCD spectra were obtained with the 
IGeoE 0.356-m SC telescope working at F/ll. The spectro- 
graph is a Littrow LHIRESIII with a 2400 grooves per mm 
grating and a spectral resolution of about 14.000. The initial 
reduction (bias subtraction, flatfielding, creation of 1-D im- 
ages, and wavelength calibration) of the spectra was made 
by JR. 

The rectification and removal of cosmics and flaws of all 
spectrograms were carried out in a uniform way by JN in 
SPEFO. The program SPEFO was also used to RV measure- 
ments, based on a comparison of direct and flipped images of the 
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Fig. A.l. A comparison of independent RV measurements of the 
steep Ha emission wings (upper panel) and shell core absorption 
(bottom panel). 

spectral line profiles. Since we were searching for small RV vari- 
ations and since the setting on the steep wings of the emission- 
line profiles was not always straightforward (see below), these 
RV measurements were carried out independently by JN and PH. 
Besides the settings on the steep wings of the Ha emission, we 
also measured the Ha absorption core on all spectra where such 
ab sorption was present t o have a comparison with the results 
of Katahira et al. (1996b). We also tried to measure RV of the 
He I 6678 A absorption wings but due to weakness of this line 
and its possible structure, these measure ments turned out to be 
useless so we did not use them. Following iHornetalJd 19961) . we 
also measured selected stronger and unblended telluric lines in 
all spectra and used them to a correction of the RV zero point. 
Thanks to that, the spectra from all observatories can be treated 
as coming from one instrument for all practical purposes. 

A comparison of the two sets of independent RV measure- 
ments is shown in Fig. IA. 1 1 In general, the agreement is good. A 
formal regression between the measurements of PH and JN was 
derived. Its slope is 0.98 + 0.01 for the emission and 0.94 + 0.01 
for the absorption. For the absorption line, it is conceivable that 
in specific cases one or the other measurer was confused by a tel- 
luric line blended with the stellar absorption core. For analysis, 
we used the mean RVs of the two independent measurements. 
All our RVs with the corresponding HJDs of their mid-exposures 
are provided in Table I A. fl 
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Table A.l. Radial velocities of the Ha emission wings and shell absorption core obtained via averaging the independent mea- 
surements by J. Nemravova and P. Harmanec; DAO = Dominion Astrophysical Observatory, Victoria; ROZ = Rozhen National 
Observatory ; OND = Ondfejov Observatory; LIS = IGeoE-Lisbon; OHP = Haute Provence Observatory; 



Time of obs. 


KV(Ho' em.) 


D\7/U^, ,. U. \ 


Source 


Time of obs. 


D\77U^., am \ 

KV(Jrla em.) 


D\7/U^, \ 


Source 


/tuti~"\ o/innnnn\ 
(HJD-24UUUUU) 


[Km s J 


[Km S J 




(HJD-24UUUUU) 


[Km s J 


[km S J 




4QS81 

tyJO 1 .JO / J 


4 04 




ONF) 

U11J7 


JZ7J / . JZ7U 


8 80 
o.ou 


-0 78 


OND 


4Q(S^4 (S781 


76 

u. / u 




ONF) 
uini7 


S7Q78 7Q(S0 

JZ7 / O.Z7UU 


86 
7.0U 


06 
u.uu 


ROZ 


4Q(S44 S4^Q 


-1 06 

J.7U 




ONF) 


S7QQ7 4704 


51 

y. j j 


5 77 

J. / z. 


OND 


4Q(S58 47^7 


-0 54 




ONF) 


S^077 ^770 

J JUZ / . J / l\J 


8 16 

O. JU 


2.20 


OND 


4Q(S"SQ 4898 


0.95 




OND 


S^07Q 7Q71 

JJUZ7.Z7 / 1 


9.82 


3.96 


OND 


4Q(S61 44 S 5 


1.21 




OND 


S^047 771 7 


8 00 

0.7U 


1 48 


ROZ 


4Q(S67 S708 


0.34 




OND 


S^047 777fS 

JJUH^Z.Z/ /U 


8 71 
0. / j 


74 

U. / T - 


ROZ 


T7U/7. JtJU 


2.32 


-10 51 

±U. J J 


OND 


JjUtt.Z7UJ 


10.57 


4 86 
7.0U 


ROZ 


4978fS 717Q 

^7 / OU. t \- 1 J 


2.17 


-0.69 


DAO 


51044 101 7 

JJU PT.JUl / 


10.35 


4.64 


ROZ 


40010 5571 

77:7 JU — ' J / J 


7 16 

/ . 1U 




U11J7 


51046 1601 


16 
y. ju 


1 00 

J.U7 


ROZ 


40048 5677 


74 




ONF) 


51048 1700 


1 1.01 


6 11 

U. J J 


OND 


40040 (SITS 7 


8 08 




ONF) 


51060 7807 


50 

7.77 


5 48 


OND 


50001 5117 

juuu i .J J J / 


9.63 




OND 


51087 7765 

JJUOZ.Z / UJ 


10.21 


8.46 


OND 


5001 s 4.5no 


8 06 




ONF) 

U11J7 


51101 7588 


65 

7.UJ 


7 68 


ROZ 


501 04 41 17 


4.68 




OND 


51716 5807 

JJZ.1U.JOU / 


7.69 


2.99 


OND 


SOI 77 ^7Q8 

JUl ZZ. JZ70 


8 05 

O.UJ 




OND 


JjZlU.JOtl 


6 77 
u. / ^. 


2.61 


OND 


501 50 1744 

JUJ » 7 . J Z 1 1 


7.52 




OND 


51716 5188 

JJZJU. J JOO 


7.07 


0.62 


OND 


■sn^ i fs 511 i 


12.67 




U1>J7 


51716 541 8 


6 67 
u.u / 


1 60 

1 .U7 


OND 


50410 5517 


12.17 




v7 - > 1 7 


51716 5444 


6 85 

U.OJ 


7 75 

z. / J 


OND 


50410 1454 

JU^t J7. J7 J7 


1 1 07 






51744 5817 


7 85 

/ .OJ 


17 

U. JZ 


ROZ 


50448 4556 


1 7 10 




ONF) 


51101 4677 


6 80 


-1 16 
J. 1 u 


ROZ 


50508 1166 

J U J UO . J J UU 


58 

7.JO 




ONF) 


51101 4714 


6 76 
u. / u 


-2.61 


ROZ 


50500 17 11 

J UJ U7 . J Z 1 J 


8.40 




OND 


51106 4678 

JJJUU.tUZO 


6.94 


-2.43 


ROZ 


S1 777 (S8^4 

J 1ZZ / 


1 5 86 




DAO 


^Q^9 70^4 

JJJJZ.Z7Jt 


6 14 
u. n - 


-1 80 

J.OU 


ROZ 


S1481 S805 


1 6 66 




v7 - > 1 7 


JJJJJ.JZ70 


4 87 


-4 51 


OND 


S1 S70 7QS4 


1 8 00 

10.U7 


-1 1.26 


OHP 


JJJJJ.J UJ 


4 1 5 
^t. 1 j 


-5 51 

j. j j 


OND 


51 577 7711 


15.23 




OHP 


S^4S7 747S 


4.60 


-5.45 


ROZ 


51 577 7871 


15.27 


-7 80 


OHP 


51555 5651 


4 75 


-6 60 
-u.uu 


OND 


51 571 1101 

J J- J /J.J1U1 


14.43 




OHP 


51570 5661 

J J J / 7.JUUJ 


4 07 


-0 1 5 

7.1 j 


OND 


51 571 11 84 

J 1 J / J. J 1 ot 


14 80 




OHP 


5161 5 5678 
jjuij.juz.0 


7 54 

z,. J^t 


-6 85 

-U.OJ 


OND 


51 576 7801 

Jl J / U. ZOU 1 


20.66 




OND 


51678 0077 

JJUZO.UU/ / 


3.98 


-4.82 


DAO 


51707 5007 

Jl / 7J — '7Ui 


14.92 


-1 58 




51618 0848 


4 01 


-5 06 


DAO 


51888 4718 


15.65 


-12.79 


OHP 


51651 5018 

J JUJ 1 .J 7 JO 


2.64 


-4.98 


OND 


51880 4170 

J lOOy.T^l ZU 


16.89 


-1.85 


OHP 


51658 5115 

-J JUJ O.J 1 J J 


3.02 


-6.07 


OND 


51807 4705 

J J. G 7i-. i / U J 


16.33 


-3.26 


OHP 


51708 8487 

J J / UO . UtOZ 


3.65 


-2.43 


DAO 


51075 1171 
j 1 7 / j. j jz± 


14.73 


-7 08 


OND 


51708 8557 

JJ / UO.OJJi 


1 68 

J.UO 


-4 75 


DAO 


57716 4148 


12.34 


-0.05 


OHP 


51745 4108 

J J /tJ.tJUO 


3.12 


-3.33 


OND 


57716 411 1 

JZZJU.tJ u 


1 1 40 


-0 60 


OHP 


S^74S 4444 
jj/ 1 j .tttt 1 


5 75 

j.zj 


-1 76 

- J.ZU 


OND 


57717 4101 


13.64 


1.47 


OHP 


^7R7 ^4S^ 

J J / OZ.JtOJ 


-6.17 


-13.71 


ROZ 


57718 4717 

JZZ JO.T^i 1 Z 


1 5 07 

1 J .u 1 


-2.16 


OHP 


51701 1470 

JJ / 7l.J t tZ.U 


-6 16 
u. 1 u 


-0 60 

7.U7 


OND 


57710 4008 


14.07 


5.44 


OHP 


51701 6055 

JJ / 71.U7JJ 


-4.84 


-10.42 


DAO 


57741 4048 


13.66 


2.18 


OHP 


5181 1 6687 

J J O 1 J .uuoz 


0.54 


-6.90 


DAO 


S7741 41 (S8 

JZ.Z.H 1 1 1 UO 


1 1 78 

1 J.ZO 


-0.09 


OHP 


JJolJ.U / Jl 


10 

U. 1U 


-7 00 

/ .uu 


DAO 


57747 1068 

JZ.ZtZ. J7UO 


13.50 


-2.44 


OHP 


51814 1110 

JJO it.JJ 17 


-2.30 


-11.32 


ROZ 


S7747 4088 


13.17 


-1.56 


OHP 


JJOlt. / UJO 


15 

U. J J 


-5 74 
- j. / 7 


DAO 


57761 1571 


4.86 


-21.12 


OHP 


J J O 1 7 .zouu 


1.72 


-3.85 


OND 


S7764 ^5Q8 

JZ.zUt.JJ70 


1 00 

j.77 


-20.28 


OHP 


JH-UUZ.UUUZ 


-5 50 

-J. JU 


-5 17 

j. 1 / 


DAO 


S7764 ^7SQ 

JZ.ZUt.J / J7 


4 10 
t. 1 y 


-20.14 


OHP 


S404Q 40^^ 

J i tU t T7 .tUJ J 


1 00 
1 .uu 


-0 40 


ROZ 


S7787 77S^ 

JZ.Z.O / . / / J J 


12.91 


-1 06 


DAO 


S404Q SASR 

J t tU t t7 . J U O O 


50 

U.J7 


-4 17 

-t t. JZ 


ROZ 


JZJJJ.UJUU 


1 7 70 

li. / U 


-0 50 

U.J7 


DAO 


JtUJ 1 .JO J / 


-0 08 
u.uo 


-0 54 


ROZ 


S7(S64 7S87 

JzUUt.zJO / 


1 7 66 


6.15 


OHP 


JtUOJ . 1 7 JU 


4 10 

T". JU 


6 64 
u.ut 


OND 


S770(S 7^7S 

JZ / UU. / J / J 


9.20 


3.07 


DAO 


S40Q7 ^W7 

J i +U7 / .JjJJjL 


5.58 


5.29 


OND 


S7710 74^4 

JZ / lU.Z.tj'T 


8 04 


-0 88 
u.00 


ROZ 


S41 O'S 7S^7 

Jtl U J . / J JZ 


6 78 


4 50 


DAO 


S7860 S7Q4 

JzOUU. JZ.7H 


1 1 01 

1 j.u 1 


5 71 

j. / 1 


ROZ 


Jtl UO.JJUt 


4 60 


1 58 

J.JO 


ROZ 


S7877 S8^S 

JZO / / . JOJJ 


14.12 


4.46 


OND 


S41 1 S ^7^0 

Jtl 1 J.JZJU 


50 

7.77 


7.83 


OND 


57800 5051 

J ZO/7 .J7 J J 


9.18 


2.73 


OND 


541 1 5 1165 

jTllJ . J JU J 


6.24 


5.07 


OND 


S7Q00 SfS77 

JZ.7UU.JU / z 


10 16 

1U. 1U 


1.06 


OND 


S41 1 ft ^0Q7 

Jtl 1U.JU7Z 


7 04 

/ .77 


5 01 

J.Ul 


OND 


S7Q00 S6Q4 

JZ7UU. JU7H 


01 

7.71 


1 46 

J.TU 


OND 


S41 17 ^7Q1 

Jt 11/ . JZ7 1 


8 40 
0.77 


4 07 


OND 


S7Q07 4^77 


1 1 08 


71 

U.ZJ 


OND 


S41 76 7777 

J^+lZU.ZZZ / 


5 71 


1 45 

J.7J 


OND 


S7Q07 4^76 

JZ7UZ.H-J / U 


10 01 

1U.U 1 


84 


ONF) 

V7 - \ 1 7 


S41 7844 

JH- 1 J J.ZOH-H- 


10 11 

1U. 1 J 


55 

7.J7 


WIN U 


52904.6492 


10.21 


-1.22 


OND 


54162.3401 


10.00 


7.07 


OND 


52904.6537 


8.10 


-1.58 


OND 


54164.2850 


9.44 


8.85 


OND 


52904.6581 


8.64 


-2.50 


OND 


54186.2921 


8.34 


10.40 


OND 


52949.6009 


4.53 


-4.85 


OND 


54188.3212 


8.24 


7.65 


OND 


52949.6044 


7.99 


-1.69 


OND 


54341.0107 


8.10 


9.45 


DAO 


52949.6091 


8.48 


0.25 


OND 


54387.5177 


7.60 


9.07 


OND 


52952.5541 


9.42 


1.66 


ROZ 


54442.9976 


-1.19 


7.27 


DAO 


52952.5619 


8.61 


0.73 


ROZ 


54490.2613 


2.53 


5.75 


OND 


52957.5042 


10.46 


1.38 


OND 


54490.9272 


5.12 


7.32 


DAO 


52957.5113 


9.64 


-0.03 


OND 


54508.3482 


5.75 


6.33 


OND 



J. Nemravova et al.: Duplicity of BU Tau, Online Material p 3 

Table A.2. (cont.) Radial velocities of the Ha emission wings and shell absorption core obtained via averaging the independent 
measurements by J. Nemravova and P. Harmanec; DAO = Dominion Astrophysical Observatory, Victoria; ROZ = Rozhen National 
Observatory ; OND = Ondfejov Observatory; LIS = IGeoE-Lisbon; OHP = Haute Provence Observatory; 



Time of obs. 


KViHa' em.) 


KV(iiQ' aus.j 


Source 


Time of obs. 


K\(Ha em.) 


KV(Hff abs.) 


Source 


(HJD-24UUUUU) 


[Km s J 


[Km s J 




(HJD-24U0U0U) 


[Km s J 


[km S J 




S4S 1 Q (S647 


7 47 


7 46 


DAO 


S4871 4.085 


1 45 


87 

\J.O 1 


ONn 

U1^17 


S4ST7 ^14S 
Jtjj / .j itj 


7 61 


7 0? 


v7 - > 1 7 


^4871 4^47 


-1 05 


-1 ^5 






7.JO 


7 5^ 


WIN l-' 


5487? 7^1 ^ 


-0 50. 


00 

U.V7U 


ONn 


S4SS7 ^17Q 


7.UU 


10 8fi 




54877 757^ 


-0 47 


-0 18 


ONn 


S471 R S1 ^8 

Jt / 1 O.J 1 JO 


Q "3,4 


8 46 




^4874 14^56 


-4 1 5 


-7 85 


T TS 

I^IO 


S474S df\]f\ 


X 77 

o. / / 


8 77 
o. / / 


OND 


54880 ^HlSO 


-1 Q1 

1.71 




LIS 


5475'? AAS7 


10.20 


9.03 


OND 


54881 3730 


-2.97 


2.77 


LIS 


S47S^ 4S40 


Q 06 


8 If, 
o. / u 


Awn 

WIN 17 


S4RS? "Kill 

JH-OOZ. JZZZ 


^ ^7 


1 51 


l^lo 


54761.4009 


10.36 


9.77 


OND 


54911.6519 


5.48 


6.32 


DAO 


54763.4827 


10.24 


8.19 


OND 


54911.6904 


5.42 


6.43 


DAO 


54798.4322 


10.35 


9.47 


OND 


54912.6793 


5.60 


6.45 


DAO 


54804.2374 


9.88 


8.46 


OND 


54924.3033 


7.18 


7.77 


OND 


54840.4814 


8.89 


9.77 


OND 


55050.5248 


11.68 


8.25 


OND 


54857.4154 


8.08 


7.49 


OND 


55071.5298 


11.44 


7.33 


OND 


54862.6448 


5.96 


3.92 


DAO 


55083.6501 


6.07 


3.72 


OND 


54862.6812 


5.78 


3.76 


DAO 


55097.4911 


-1.45 


4.41 


OND 


54863.6276 


5.43 


4.58 


DAO 


55112.4081 


3.52 


-3.81 


OND 


54863.6627 


4.72 


4.21 


DAO 











